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FRODGE.  REMONDI  &  LAPUCHA 

Real-Time  Centimeter  Positioning  with  GPS: 

A  System  Developed  by  the  US  Army  Corps  of  Engineers 

Ms.  Sally  L.  Frodge,  US  Army  Topographic  Engineering  Center 
Dr.  Benjamin  Remondi.  National  Geodetic  Survey 
Or.  Dariusz  Lapucha,  John  E.  Chance  Associates,  Inc. 

INTRODUCTION 

The  U.S.  Army  Corps  of  Engineers  (USACE)  is  responsible  for  keeping  the  waterways  of  the 
United  States  of  America  navigable.  Recognizing  the  need  for  and  potential  benefit  of  a  more  accurate 
positioning  system  for  its  hydrographic  surveying  and  dredging  mission,  the  Corps  embarked  on  an 
ambitious  research  program  to  develop  a  prototype  positioning  system  based  on  Global  Positioning 
System  (GPS)  carrier  signals.  This  prototype  system  is  designed  to  deliver,  in  real  time,  three- 
dimensional  positions  with  subdecimeter  accuracies,  for  ranges  up  to  20  kilometers  (km)  using  a  single 
reference  station.  This  project  is  in  its  final  year  of  a  planned  6  year  effort  that  has  gone  through  a 
full  development  cycle  from  concept  to  development  of  a  working  prototype.  Feasibility  studies, 
system  analysis,  extensive  testing  strategies,  and  considerable  research  with  regard  to  the  resolution 
of  carrier  ambiguities  'on-the-fly'  were  some  of  the  major  program  activities. 

Several  millions  of  dollars  are  spent  each  year  by  the  Corps  in  the  maintenance  of  the 
waterways  of  the  United  States  of  America.  The  Corps  conducts  Condition  Surveys  on  a  routine  basis 
to  identify  channel  obstructions.  If  dredging  is  required  to  clear  the  channel  obstruction,  the  Corps 
often  contracts  this  work  out.  If  this  is  the  case,  a  typical  scenario  is  for  the  Corps  to  perform  a 
hydrographic  survey  before  the  work  is  done  to  identify  to  the  contractor  the  specific  material  to  be 
removed.  After  the  work  is  completed,  the  Corps  performs  a  Contract  Payment  survey.  The  current 
horizontal  accuracy  standard  used  for  a  Contract  Payment  survey  is  3  meters  (m)  1  ORMS.  The 
majority  of  the  positioning  systems  used  today  by  the  Corps  are  range-range  or  range-azimuth  systems. 
Most  of  these  systems  require  daily  calibration  or  initialization  at  a  site  local  to  the  job.  In  addition, 
occupation  of  other  previously  surveyed  shore  stations  are  necessary  for  each  day  of  the  survey. 

Survey  and  dredging  operations  also  must  be  referenced  to  a  vertical  datum.  Tidal,  river  or  lake 
gages  are  used  to  establish  readings  upon  which  a  vertical  datum  such  as  mean  lower  low  water 
(MLLW)  is  established.  The  surveyor  must  determine  the  differences  in  elevation  from  the  established 
datum  using  methods  such  as  zoning  models,  which  typically  limit  the  accuracy  to  0.2  m.  Rnal 
payment  is  made  to  the  contractor  when  both  the  contractor  and  the  Corps  agree  that  the  identified 
material  was  removed  and  on  the  volume  of  the  removed  material.  The  prototype  system  being 
developed  by  the  Corps  is  designed  to  provide  subdecimeter  accuracies  in  three  dimensions  in  real¬ 
time.  This  system  does  not  require  static  initialization:  initialization  is  on-the-fly  (OTF)  and  typically 
within  1 5  seconds  or  less.  Only  a  single  shore  station  is  necessary;  daily  calibration  is  not  necessary. 
It  is  anticipated  that  the  implementation  of  this  system  throughout  the  Corps  will  save  the  government 
and  its  taxpayers  millions  of  dollars. 

BACKGROUND 

The  development  of  the  On-The-Fly  (OTR  prototype  system  was  initiated  by  the  U.S.  Army 
Topographic  Engineering  Center  (TEC  )  in  1988  and  funded  by  the  Dredging  Research  Program.  The 
project  progressed  in  three  primary  phases  which  have  been  well  documented  in  previous  papers.^-’'* 
Dr.  Remondi  of  the  National  Oceanic  and  Atmospheric  Administration  /  National  Geodetic  Survey 
(NOAA/NGS)  developed  the  base  algorithms  for  the  OTF  system  developed  through  phase  two.^‘*‘*''° 
The  goal  of  phase  three  was  to  develop  a  working  prototype  that  would  be  ready  for  demonstration 
no  later  than  October  of  1 993.  In  order  to  accomplish  this  goal  and  meet  the  schedule,  a  contract  was 
signed  with  John  E.  Chance  Associates  (JECA)  to  achieve  an  operational,  real-time  prototype  system 
from  the  current  post-processing  research  and  development  capability.  A  working  version  of  OTF 
post-processing  software  was  developed  by  Dr.  Remondi  during  phase  two  and  was  used  to  analyze 
the  data.  The  third  generation  of  this  software  was  completed  by  October,  1 992,  and  the  modification 
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of  this  version  eventually  became  the  real-time  OTF  software  prototype.  The  movement  of  the  post¬ 
processing  versions  onto  a  real-time  platform  and  interfacing  the  hardware  components  into  a  working 
prototype  was  a  significant  challenge.  Static  baseline  testing  of  the  initial  prototype  system  began  in 
March  of  1 993  at  the  JECA  facility  in  Lafayette.  Louisiana.  Several  iterations  of  this  system  were 
accomplished  and  have  been  tested  and  improved  since  the  spring.  Phase  three  has  been  completed 
and  the  program  has  moved  into  the  demonstration  phase.  Except  for  the  very  final  stages  of  /9eta 
testing,  the  project  has  progressed  to  completion.  The  first  demonstration  of  the  OTF  prototype  was 
in  August,  1 993,  with  the  first  public  demonstration  following  in  October.  Demonstrations  have  taken 
place  at  several  locations  throughout  the  USA  (Virginia.  North  Carolina,  Oregon,  California  and  Texas). 
International  demonstrations  have  included  Scotland.  Australia  and  England.  The  OTF  system  has  been 
used  as  a  production  level  engineering  system  on  several  jobs  throughout  the  USA  (Montana,  New 
Mexico,  Texas,  and  California).  Most  recently,  the  OTF  system  was  installed  and  tested  onboard  a 
sea  going  dredge.  Performance  of  the  OTF  system  has  exceeded  expectations. 

THEORETICAL  BASIS  FOR  KINEMATIC  ON-THE-FLY  GPS 

The  development  of  the  algorithms  necessary  for  high-precision  kinematic  On-The-Fly  GPS 
(KOTF  GPS)  began  in  late  1989  and  is  documented  in  numerous  publications. Initialization  is 
performed  in  three  steps:  First,  a  meter-level,  first-guess  position  determination  is  achieved;  second, 
a  search  grid  is  established;  third,  the  grid  "candidates*  are  evaluated  to  isolate  the  correct  grid  point. 
In  fact,  there  are  many  variants  and  the  following  approach  is  both  representative  and  simple  to 
describe. 


Stan  1 :  Meter  Level  Positional  Boundary 

The  initial  estimate  of  the  position  of  the  user's  receiver  can  be  made  from  the  code  range 
measurements.  The  equations  are: 


Where; 
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is  the  code  range  measurement  u 

is  the  range  model  in  meters  b 

is  the  unmodeled  part  n 

represents  the  other  satellites 
needed  to  form  the  double  difference  k 
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:  represents  the  user  receiver 
:  is  the  base  or  reference  receiver 
:  represents  one  less  than  the  total 
number  of  satellites 
:  represents  the  reference  satellite 


It  is  assumed  that  there  are  at  least  three  satellites  other  than  k;  most  of  the  time  there  will  be  four 
to  six  other  satellites  (i.e.,  j  =  j,,  y'j, ...  ,i^  where  n  =  4,  5,  or  6).  Although  this  system  of  equations 
can  be  solved  for  the  position  of  user  receiver  (i.e.  p  =  p  (  zj  better  results  are  achieved  by 
smoothing  the  code  ranges  with  carrier  range  measurements.  The  carrier  range  equation  is: 

A  ( +  N^*^  ]  =  pj-  +  e. 


Where: 

A  :  is  the  wavelength  of  the  carrier  signal  (meters  per  cycle) 

ip  :  is  the  carrier  phase  measurement  (cycles) 

N  :  is  the  unknown  carrier  phase  integer  ambiguity  (cycles) 

:  is  the  unmodeled  part  of  the  carrier  (meters) 

The  same  range  model  p  is  in  both  equations.  Since  N  is  a  constant  bias,  N  will  drop  out  of  a  time 
difference,  leaving  the  following: 

R'",^,  (to)  =  it,)  -  A[  Vbu*  It.)  -  <PJ^  (to)  1  + 

Simply  put,  code  range  measurements  at  any  subsequent  epoch,  t,,  can  be  mapped  to  the  reference 
time  (epoch),  („  .  This  provides  a  large  number  of  different  measurements  of  (tj  which  can  be 
averaged.  Finally,  placing  these  averaged  ranges  into  the  initial  code  range  equations  allows  the  ( 
zj  to  be  determined  at  the  meter  level.  The  search  grid  of  step  two  is  built  around  this  initial 
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positional  estimate.  The  better  the  initial  positional  estimate,  the  smaller  the  search  volume  can  be. 
This  equates  to  a  shorter  initialization  time  since  many  otherwise  attractive  grid  candidates  will  be 
eliminated  due  to  the  fact  that  they  fall  outside  the  boundaries  of  the  search  volume.  Stated  another 
way,  the  smaller  search  volume  leads  to  faster  computational  times  since  fewer  candidates  need  to 
be  considered. 

Step  2:  Forming  the  Search  Grid 

There  are  many  ways  to  form  the  search  grid.  One  efficient  way  is  based  on  the  intersection  of 
three  double  difference  planes  from  a  given  set  of  four  satellites.  Dr.  Remondi  has  implemented  two 
other  search  schemes  (not  presented  here),  which  do  not  depend  on  the  intersection  of  three  double 
difference  planes.  By  selecting  four  satellites,  each  of  which  is  not  too  low  in  elevation  and  which 
together  provide  a  favorable  Positional  Dilution  of  Precision  (PDOP),  the  real  values  of  carrier  phase 
ambiguities  can  be  computed  at  the  initial  positional  estimate  and  rounded  to  closest  integers.  Placing 
these  determined  integers  into  the  carrier  phase  equations  for  those  four  grid  satellites  results  in  a 
single  grid  point.  Neighboring  grid  points  can  be  computed  by  incrementing  any  one  of  the  three 
integer  ambiguities  by  unity.  Having  found  three  orthogonal  neighboring  grid  points,  others  can  be 
computed  within  predefined  neighboring  volumes  by  vector  addition.  Significant  efficiencies  can  be 
achieved  by  performing  this  procedure  for  both  L,  and  L,,  although  in  principle  one  can  do  it  for  just 
L,  (i.e.,  while  still  using  L,  and  L,  data  in  the  evaluation  phase.)  The  points  of  each  grid  are  formed 
as  the  intersection  of  three  double  difference  planes  and  thus  the  grids  are  actually  a  lattice  of  three 
dimensional  positions.  There  will  be  a  limited  number  of  grid  points  that  will  appear  in  close 
proximity  to  an  L,  grid  point  identified  through  the  process  outlined  above.  These  remaining  grid 
candidates  can  be  evaluated  to  determine  which  one  is  the  correct  candidate.  For  example,  if  there 
are  4000  L,  grid  points  and  2000  L}  grid  points  for  a  total  of  6000  grid  points,  somewhere  around  200 
candidates  will  remain  to  be  processed  in  the  evaluation  step  described  below. 

Step  3.  Evaluating  the  Candidates 

This  is  the  final  step  of  the  KOTF  initialization  process.  For  a  test  grid  candidate,  the  carrier 
phase  equations  are  used  to  determine  the  integer  ambiguities  for  all  double  differences.  This  permits 
a  computation  of  the  modeled  range,  p,  and  ultimately  provides  a  residual.  The  correct  grid  candidate 
will  have  small  residuals  whereas  the  others  will  not.  Should  multiple  grid  candidates  have  similar  small 
residuals,  lane  resolution  may  not  be  possible.^  The  calculated  statistics  will  indicate  clearly  if  a  correct 
candidate  has  been  found,  and  that  the  correct  solution  for  the  integers  has  been  resolved.  If  the 
resolution  is  uncertain  at  the  reference  time,  tg,  one  can  continue  the  process  in  subsequent  epochs, 
tf.  until  a  single  lane  clearly  emerges  as  the  correct  candidate. 

Great  effort  has  been  expended  towards  determining  reliable  statistical  acceptance  criteria  for 
the  OTF  ambiguity  resolution  procedure.”-'  Stated  simply,  how  much  confidence  does  one  have  that 
correct  lanes  have  been  isolated.  The  chosen  statistical  norm  is  based  on  a  mean  of  absolute  values 
of  the  computed  real-valued  double  difference  ambiguities  from  their  rounded  integer  values.^  The  OTF 
ambiguity  resolution  is  accepted  if  it  gives  at  least  3  percent  discrimination  between  the  top 
candidates.  This  statistical  criterion  yields  roughly  99  percent  success  rate  of  the  real-time  system 
ambiguity  resolution.  When  a  still  higher  level  of  certainty  is  required  5  percent  is  used. 

DESCRIPTION  OF  THE_QTF  SYSTEM 

The  real-time  OTF  prototype  delivers  high  precision  kinematic  positioning  accurate  to  less  than 
5  cm  at  the  antenna  phase  center,  while  simultaneously  providing  a  separate  output  for  meter  level 
differential  GPS  (DGPS)  for  navigation  purposes.  The  heart  of  the  system,  the  OTF  software,  was  the 
focus  of  much  of  the  development  effort.  The  software  was  designed  to  require  minimum  operator 
attention  and  has  several  quality  control  procedures  built-in  to  assure  that  the  highest  performance  and 
reliability  of  the  system  is  maintained.  The  system  consists  of  setups  at  the  remote  and  reference 
stations  and  a  data  link  as  shown  schematically  in  Rgure  1 .  Dual  frequency  (Lf/Lj)  GPS  receivers  are 
required  at  both  stations.  Although  the  primary  goal  was  the  real-time  capability,  post-processing 
software  was  also  developed  that  serves  as  an  excellent  engineering  tool. 

Shipping  the  raw  data  from  the  reference  station  allows  a  single  386  SX  Personal  Computer  (PC) 
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to  carry  out  the  reference  station  functions,  if  the  original  configuration  as  depicted  in  Figure  1  is  used. 
These  functions  include  setting  the  GPS  receiver  to  output  the  required  data,  translating  that  data  to 
the  desired  format,  and  transferring  the 
formatted  data  to  the  data  link  for  transmission. 

Additionally,  the  reference  station  package  is 
capable  of  recording  raw  GPS  measurements  if 
the  operator  requests  it.  Although  the  prototype 
has  been  developed  and  built  using  Intel-based 
386  and  486  computers  and  Trimble  4000  SSE 
receivers,  it  is  hoped  that  in  the  near  future  other 
platforms  and  receiver  types  may  be  interfaced 
to  use  the  OTF  software.  The  system  has  been 
modified  from  the  original  configuration  to  allow 
the  user  to  select  the  option  of  moving  the 
processing  functionality  to  the  monitor  station 
site.  This  latter  configuration  requires  486  DX 
computers  at  both  station  sites  but  allows  the 
user  to  collect  data  at  the  remote  site  and 
process  the  data  at  the  base  station.  This 
configuration  was  implemented  for  an  upcoming  Figure  1.  Block  Diagram  of  the  OTF  Prototype 
axperiment  that  will  track  a  floating  buoy  that  System's  (Original  Configuration) 
will  be  collecting  data  to  pursue  real  time 
modeling  of  tides.  Since  the  buoy  cannot  be 

occupied  as  a  vessel  can,  the  data  are  sent  back  to  the  base  reference  site  for  processing  onshore. 

High-precision  KOTF  positioning  is  available  from  the  system  once  integer  ambiguities  are 
resolved  by  the  software:  five  satellites  are  required  for  this  process.  As  long  as  the  system  remains 
in  the  KOTF  mode,  real-time  subdecimeter  positioning  in  three  dimensions  is  available  at  the  remote 
site.  Remaining  in  the  KOTF  mode  requires  both  reference  and  remote  station  data  as  well  as  a 
maintenance  of  lock  on  at  least  four  satellites.  If  that  number  drops  to  below  four,  the  ambiguities  will 
again  be  resolved  after  the  system  reacquires  lock  on  five  satellites.  The  software  is  "smart”  in  that 
it  will  automatically  detect  the  need  to  reinitialize.  The  software  also  will  trigger  reinitialization  if 
quality  factors  based  upon  residuals  fail  to  meet  certain  predefined  limits.  Note  that  the  system  is 
always  capable  of  meter-level  DGPS  navigation  even  when  loss  of  lock  occurs.  Furthermore,  the 
system  will  provide  this  navigation  function,  for  a  limited  time,  even  without  data  from  the  reference 
station. 

The  system  uses  Lf/Lj  carrier  phase  and  C/A  code  ranges  for  ambiguity  resolution.  The  system 
has  been  designed  not  to  rely  on  the  continuity  of  L,  carrier  phase,  since  only  the  fractional  phase  part 
of  the  Lj  carrier  measurements  is  used  in  the  OTF  ambiguity  resolution.  After  ambiguities  are  resolved, 
only  the  L,  carrier  ranges  are  required  to  maintain  the  high-precision  KOTF  positioning.  The  meter-level 
OGPS  process  uses  primarily  C/A  code  and  L,  carrier  ranges.  Within  the  OTF  software,  the  time 
interval  for  the  initialization  process  is  user  selectable.  Typically,  this  time  is  set  at  15  seconds,  (15 
1  -second  epochs)  although  several  of  the  initial  tests  over  the  summer  were  run  using  a  30  second 
time  interval.  The  data  are  first  analyzed  in  a  forward  manner,  from  the  reference  epoch  forward  in 
time  towards  epoch  t„.  if  the  OTF  system's  software  cannot  resolve  the  integers,  then  the  software 
automatically  will  begin  the  initialization  process  again,  processing  the  same  data  backwards  in  time, 
beginning  with  epoch  t„  and  proceeding  backwards  towards  t„.  Although  the  initial  1 5  second  forward 
pass  is  normally  sufficient,  in  situations  when  additional  time  segments  are  required,  initialization  may 
not  require  the  entire  set  of  epochs  back  to  to,  i.e.,  initialization  may  take  more  than  1 5  seconds,  but 
less  than  30.  The  software  will  indicate  if  initialization  was  obtained  on  the  forward  or  backward  pass. 
The  computational  time  required  for  OTF  ambiguity  resolution  is  usually  less  than  one  second. 

The  prototype  (in  its  original  configuration)  was  developed  on  the  premise  that  all  required  raw 
GPS  observations,  i.e.,  the  GPS  time  tag,  Lf/Lj  carrier  phase  and  L,  code,  are  transmitted  from  the 
reference  station  and  the  actual  computations  necessary  for  KOTF  and  code  OGPS  take  place  at  the 
remote  processing  site.  The  KOTF  process  requires  time-matched  reference  and  remote  station  data. 


In  the  code  differential  process,  extrapolated  reference  station  differential  corrections  are  applied  to 
the  current  remote  station  observations,  as  is  done  in  standard  OGPS  systems.  The  difference  here 
is  that  these  corrections  along  with  their  rates  are  generated  at  the  remote  site  from  the  raw  data  that 
was  received  from  the  reference  station.  As  mentioned,  the  processing  function  can  now  be  carried 
out  at  either  the  reference  or  remote  station.  The  system  can  be  interfaced  with  any  other  system 
requiring  three  dimensional  centimeter  accuracy  in  positioning  using  the  output  interface  string 
developed  for  this  project.  This  string  was  designed  to  be  as  close  as  possible  to  current  National 
Marine  Electronics  Association  (NMEAI  0183  GPS  string  formats,  a  widely  accepted  standard  for 
marine  applications. 

The  system  requires  a  data  link  capable  of  at  least  of  4800  baud.  The  system  has  been  tested 
using  9600  Baud  UHF  sets  (460  MHz)  that  can  operate  at  line-of-sight  distances  up  to  25  km, 
depending  upon  antenna  height  and  power  used.  Tests  have  also  been  run  using  one  watt  spread 
spectrum  sets  (906-920  MHz).  Low  power  spread  spectrum  sets  have  the  advantage  of  not  requiring 
a  Federal  Communication  Commission  (FCC)  license  within  the  USA.  They  provide,  however,  a  more 
limited  range  of  possibly  10  km  (line-of-sight,  rather  ideal  conditions).  The  effective  range  obtained 
within  these  tests,  without  using  any  repeaters,  was  8.1 -8.5  km.  VHF  sets  (163.0  MHz  and  164.5 
MHz)  were  used  with  the  system  and  provided  a  disappointing  8-9  km,  but  have  since  been  retuned. 
Ranges  of  approximately  17  km  were  obtained  in  subsequent  tests  with  these  same  VHF  radios. 
Successful  tests  were  run  out  to  18.6  km  using  the  UHF  sets. 

Unlike  standard  differential  position  fixing  (i.e.,  where  the  range  and  range  rate  corrections  are 
used  to  create  a  position  fix  with  little  latency),  KOTF,  as  mentioned  before,  requires  time  matched 
data.  Any  latency  in  the  data  link  results  in  a  corresponding  delay  in  the  precise  position  computation. 
Position  latency  on  the  order  of  two  seconds  or  more  can  be  seen  on  the  KOTF  position  fix,  depending 
on  the  radio  link  used.  At  the  remote  station,  the  reference  and  remote  station  data  are  combined  to 
accurately  determine  the  remote  station  position.  The  remote  station  software,  in  its  typical  mode  of 
operation,  provides  navigation  and  KOTF  output  at  a  1  Hz  rate.  The  OTF  system  uses  off-the-shelf 
equipment  and  customized  software.  The  system  works  in  a  robust  and  reliable  manner  out  to  the 
design  goal  range  of  20  km  in  real-time. 

RESULTS  FROM  TESTING  THE  SYSTEM 

Four  primary  goals  were  targeted  for  the  testing  of  the  real-time  OTF  system:  to  determine  how 
quickly  the  correct  integer  ambiguity  for  each  satellite  can  be  established;  to  determine  the  accuracy 
of  the  real-time  30  position  measurement:  to  verify  the  repeatability  of  the  real-time  position  even 
after  reinitializing  the  satellite  integers;  and,  to  verify  that  the  system  can  determine  when  it  is  unable 
to  output  the  correct  kinematic  position.  To  meet  these  objectives,  the  testing  was  broken  into  stages 
that  progressed  from  static  baseline  tests  to  land  mobile  tests  and,  finally,  to  tests  in  an  operational 
environment  aboard  a  survey  vessel.  In  the  static  tests,  the  position  determined  in  real-time  OTF  mode 
was  compared  to  the  known  antenna  position.  In  the  mobile  tests,  a  survey  truck  was  driven  to 
existing  control  points  where  the  real-time  position  was  compared  to  known  surveyed  position.  The 
truck  would  revisit  the  points  repeatedly  with  the  system,  while  operators  put  the  system  through 
various  tests,  e.g.  instructing  the  system  to  recompute  the  satellite  integers. 

The  real-time  system  was  configured  the  same  way  for  all  of  the  tests,  whether  static,  mobile 
or  in  the  operational  marine  environment.  For  each  of  the  tests,  the  OTF  software  was  setup  as 
follows:  initialization  time  set  to  a  count  of  currently  1 5  (or  30)  one-second  epochs;  cylindrical  search 
volume  of  1 .5  m  horizontal  and  2.0  m  vertical:  L^/Lj  noise  set  to  50  mm.  The  receiver  was  set  to  track 
L,  C/A  code  and  Lj  P  code.  The  Lj  P  code  was  used  only  in  the  OTF  initialization  process  at  this  time. 
Testing  in  Lj  cross  correlation  mode  was  also  done  with  no  obvious  reduction  in  accuracy  results, 
notwithstanding  the  lower  signal  to  noise  values  which  may  require  longer  initialization.  At  both  the 
reference  and  remote  sites,  survey  grade  GPS  antennas  with  ground  planes  were  used.  During  the 
demonstrations,  there  were  some  cases  where  obstructions  coupled  with  poorer  GPS  constellation 
geometry  caused  initialization  to  require  more  than  two  passes.  Note  that  for  all  tests  and 
demonstrations,  the  system  was  turned  on  while  en  route  to  assure  initialization  occurred  under  mobile, 
dynamic  conditions. 

The  static  tests  began  in  March  of  1993.  Many  tests  were  performed  over  a  140  m  baseline 


MEASURED 

STANDARD 

MAXIMUM 

DEVIATION 

VARIANCE 

LATITUDE 

0.004  m 

0.020  m 

LONQITUDE 

0.003  m 

0.010  m 

VERTICAL 

0.010  m 

0.050  m 

Table  1  Summary  of  static  test  results. 


where  the  real-time  system  was  running 
continuously  with  the  30  real-time  position 
output  compared  with  the  known  location  of  the 
antenna,  previously  established  using  static  GPS 
survey  methods.  The  reference  station 
composed  of  a  dedicated  GPS  receiver,  laptop, 
and  radio  link  and  was  left  running  unattended. 
This  same  reference  station  was  also  used  for 
many  dynamic  tests.  The  results  from  a  typical 
static  test  are  shown  in  Table  1.  The  results 
from  the  static  test  shown  in  Table  1  are  typical 
of  what  was  seen  from  many  other  tests  run 
over  the  same  baseline.  The  statistics  shown  are  calculated  using  a  typical  24  hour  period  with  the 
system  running  through  the  full  GPS  constellation  over  that  time.  Typically,  positions  were  obtained 
98  percent  of  the  time  (23.5  hrs/day).  The  remaining  2  percent  of  the  time  integers  could  not  be 
established  due  to  poor  satellite  constellations  containing  several  low  elevation  satellites  or  not  enough 
satellites.  Results  from  monitoring  the  system  over  the  static  baseline  for  24-hour  periods,  for  months, 
show  that  the  system  delivers  three-dimensional  position  with  centimeter  accuracy. 

The  usual  setup  for  land  mobile  testing  or  demonstrations  entailed  setting  up  a  truck  or  cart  as 
the  mobile  remote  user  and  then  navigating  the  system  to  positions  previously  established  using  static 
GPS  methods.  Several  truck  tests  of  this  type  have  been  run  over  ranges,  from  the  reference  station 
to  the  mobile  remote,  varying  from  less  than  1 .0  km  up  to  1 9.5  km.  Both  short  and  long  range  tests 
produced  similar  results  for  accuracy  and  repeatability.  This  type  of  test  involved  some  stationary 
occupations  of  the  point;  these  times  were  kept  to  a  minimum.  The  system  was  also  closely  monitored 
while  en  route  between  stations.  Driving  speeds  varied  from  8-40  kph  (5-25  mph)  and  system 
performance  was  found  to  be  very  satisfactory.  Rgure  2  shows  some  typical  results  from  a  truck  test 
run  near  the  JECA  facility  in  Louisiana.  This  figure  demonstrates  that  the  horizontal  positions  check 
within  1  -2  cm  and  the  vertical  positions  within  1  -3  cm. 

The  next  step  was  to  move  the  system  onto  a  survey  vessel.  Initial  tests  were  run  in  the  vicinity 
of  Norfolk,  Virginia,  in  August,  1993.  For  these  tests,  the  OTF  reference  station  was  installed  atop 
the  USAGE  Norfolk  District  office  building  and  the  remote/user  equipment  was  installed  on  the  Survey 
Vessel  (S.V.)  Adams.  The  first  test  performed  was  the  radio  range  and  integer  initialization  test.  The 
radios  that  were  used  for  this  test  were  UHF  460  MHz  2  watt  units  running  at  9600  baud  (modem 
internal  to  radio).  The  antenna  on  the  reference  station  was  a  broadband  yagi  (directional  antenna), 
pointing  down  the  waterway  towards  the  Adams.  On  the  Adams  was  an  omni-directional  antenna, 
forward  of  the  main  mast.  The  maximum  range  was  defined  as  the  usable  distance  at  which  the  radio 
link  lost  data  for  more  than  5  seconds.  The 
range  obtained  was  18.6  km.  There  was 
no  problem  initializing  integers  at  this 
distance.  One  of  the  most  important  tests 
compared  the  real-time  vertical  positions 
obtained  over  time  using  the  OTF  system 
with  tide  gage  readings  to  determine 
vertical  accuracy  relative  to  tidal  movement 
of  the  vessel.  Additionally,  recordings  were 
made  using  a  spirit  level  as  another 
independent  check.  For  these  tests  the 
Adams  was  tied  to  the  dock  for  a  full  tidal 
cycle  of  8  hours.  Summary  results  are 
shown  in  Table  2.  OTF  positional  data 
were  logged  on  disk  and  also  recorded  with 
the  tide  gage  data,  measuring  to  a  point  on 


the  Adams.^  For  the  venical  test,  the  OTF  pjgurg  2  Kinematic  Land  Test  Results, 
reference  station  was  installed  in  a  TEC 


survey  truck  which  was  then  driven  to  one  of  the  Chesapeake  Bay  Bridge  Tunnel  islands.  These  survey 
points  had  been  previously  established  on  that  island  using  static  GPS  methods.  Two  tests  were  run: 
one  with  the  Adams  was  anchored  1 900  m  away,  and  the  other  with  it  anchored  200  m  away.  Table 
2  shows  the  results  comparing  readings  from  a  spirit  level,  tide  gage  readings,  and  the  vertical  readings 
from  the  OTF  system.  For  these  tests,  the  Adams  was  tied  to  the  dock  for  a  full  tidal  swing  of  8 
hours.  The  purpose  of  the  test  was  to  check  long-term  vertical  accuracy  relative  to  tidal  movement 
of  the  vessel.  Figure  3  shows  a  detailed  plot  of  the  data  from  the  test.  Field  notes  show  that  the 
passing  boat  traffic  coincides  with  the  larger  variations  that  can  be  seen  on  the  graph,  for  example  at 
approximately  220  minutes  and  390  minutes. 

Several  tests/demonstrations  of  the  OTF  system  have  taken  place  since  September.  The  first 
demonstration  was  the  Norfolk  test  previously  described,  which  took  place  in  August  1 993.  The  first 
public  demonstration  took  place  in  early  October  1993,  in  Wilmington,  North  Carolina.  This 
demonstration  took  place  in  the  downtown  area  just  south  of  Wilmington  on  the  Cape  Fear  River.  The 
site  was  selected  to  maximize  performance  of  the  range-range  system  that  is  normally  used  onboard 
the  Wilmington  District's  S.V.  Gillette.  For  the  demonstration,  the  Gillette  ran  several  longitudinal  and 
cross  section  lines.  Running  these  lines  took  the  vessel  underneath  a  bridge  twice,  once  on  the  way 
down  to  the  survey  area  and  once  again  on  the  return  trip.  The  demonstration  onboard  the  Gillette 
also  showed  how  the  OTF  system  performed  in  an  operational  environment,  for  example,  automatically 
reinitializing  after  experiencing  obstructions,  notably  the  Wilmington  bridge.  As  mentioned  before,  the 
system  usually  initialized  within  1 5  seconds,  although  there  were  occasions  where  multiple  initialization 
time  intervals  were  required.  Still,  collectively,  initialization  took  less  than  1  minute  after  transit 
underneath  the  bridge.  To  demonstrate  the  specific  capabilities  of  the  OTF  system,  a  cart  test  similar 
to  the  previously  described  truck  test  was  run  onshore.^ 

The  demonstration  in  Astoria,  Oregon,  took  place  in  mid-November  following  the  general  format 
of  the  one  in  Wilmington.  Software  changes,  such  as  adding  the  forward  and  backward  processing 
method  described  earlier,  were  made  between  these  two  demonstrations  to  improve  the  initialization 
time.  The  OTF  system  was  setup  on  the  Portland  District's  S.V.  Hickson.  Again,  the  demonstration 
site  was  selected  such  that  the  vessel  would  have  to  transit  underneath  a  bridge  twice  during  each 
demonstration.  No  significant  problems  were  experienced  with  the  OTF  system.  The  worst  case  for 
reinitialization  was  that  the  OTF  system  did  not  re-initialize  until  within  the  fourth  pass  (less  than  one 
minute).  The  typical  case  was  that  one  1 5  second 
pass  was  required,  and  the  system  initialized  in  the 
forward  manner  on  that  first  pass.  There  were, 
notably,  several  occasions  where  the  Hickson 
passed  under  the  bridge  and  the  system  maintained 
uninterrupted  operation  at  both  the  meter  and 
subdecimeter  level.  This  situation  did  not  occur  in 
Wilmington,  probably  since  the  bridge  at  Wilmington 
was  lower  and  some  improvements  to  the  software 
between  the  Wilmington  demonstration  and  the 
Astoria  demonstration  increased  the  robustness  of 
the  system. 

The  OTF  system  was  also  demonstrated  at  the 
1994  National  Technical  Meeting  of  the  Institute  of 
Navigation  (ION)  which  was  held  in  late  January  in 
San  Diego.  Data  collected  during  that 
demonstration  showed  similar  results  to  what  has 
been  previously  reported.  During  the  writing  of  this 
paper,  the  OTF  system  was  operating  on  the  Corps  dredge  Essayons  to  assess  system  performance 
in  an  operational  production  environment.  This  job  the  Essayons  is  working  is  10-1 5  km  off  the  coast 
of  California.  Data  are  not  yet  available  but  will  be  analyzed  for  a  variety  of  performance  factors. 
Further  testing  to  determine  vessel  dynamics  is  also  currently  taking  place. 

A  limiting  factor  for  use  of  this  system  is  the  range  of  the  data  link.  Experiments  exploring  better 
communication  links,  e.g.  satellite,  are  ongoing.  TEC  is  an  experimenter  with  the  National  Aeronautics 
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Figure  3  Tide  Gage  (  □  ):  GPS  (-f-). 
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and  Space  Administration  (NASA)  on  NASA's 
Advanced  Communications  Technology  Satellite 
(ACTS).'  The  first  phase  of  TEC's  experiment 
with  the  ACTS  was  run  in  December,  1 993,  over 
a  320  km  baseline.  Carrier  phase  data  were 
collected  in  addition  to  other  data,  but  this 
portion  of  the  data  has  not  yet  been  analyzed. 

Meter  level  DGPS  ran  over  the  ACTS  link  with  no 
problems. 

The  OTF  system  has  been  used  on 
production  level  jobs  by  JECA.  Data  gathered 
during  a  breakwater  survey  were  analyzed  to 
compare  the  vertical  component  as  measured  by 
the  OTF  system  with  that  as  measured  by  a  TSS 
heave  compensator.  The  results  are  shown  in 
Rgure  4.  This  data  sat  was  gathered  where 
there  was  a  good  deal  of  chop  to  the  water  surface.  Since  it  is  usual  that  the  vertical  errors  are  twice 
as  large  ??  >  horizontal  errors  with  GPS,  this  data  set  is  very  encouraging.  The  system  has  been 

used  on  large  scale  seismic  projects  in  Montana  and  New  Mexico  where  8,000  to  1 4,000  points  were 
laid  out  and  accurately  measured  using  the  system.  The  points  were  set  up  in  real  time  using  the 
OGPS  technique,  while  data  was  simultaneously  gathered  for  OTF  kinematic  post-processing.  Each 
surveyor  was  able  to  stake  out  approximately  1 50  points  per  day,  accurate  to  the  centimeter  level. 
This  is  significantly  faster  than  previous  methods  employed  and  additionally  provides  immediate 
verification  of  accuracy.  These  factors  equate  directly  into  cost  savings. 

PROBLEMS  FACED  DEVELOPING  THE  SYSTEM 

Design  specifications  for  the  development  of  the  prototype  dictated  that  only  off-the-shelf 
equipment  be  used  which  imposed  some  limitations  on  the  design.  A  considerable  problem  faced  by 
this  project  was  being  able  to  do  the  integer  search  quickly  enough  to  be  able  to  use  off-the-shelf  PC 
class  hardware.  In  order  to  have  a  good  starting  point  for  the  integer  search,  low  noise  L,  C/A  code 
ranges  (e.g.,  0.25  m)  are  highly  desirable  to  do  the  basic  differential  position  computation.  Advances 
in  GPS  receiver  technology  certainly  contributed  to  the  success  of  this  project.  The  ability  of  the  GPS 
receiver  to  deliver  full  wavelength  L,  carrier  phase  in  the  presence  of  Anti-Spoofing  (A-S)  provides  an 
important  processing  advantage  over  L;  squared  data.  Those  two  factors  allowed  the  project  to  be 
successful  using  a  486  DX  33  platform  for  the  remote  station  and  a  386  SX  platform  for  the  reference 
station.  Another  problem  faced  by  this  project  (and  DGPS  in  general)  is  the  current  radio  link 

equipment  that  is  available  and  which  can  be  licensed.  Both  the  capability  and  licensing  restrictions 
limit  the  extent  to  which  this  system  can  be  tested  and  used.  Currently,  about  300  bytes  per  epoch 
are  necessary  to  ship  the  L^/Lj  data  from  the  reference  to  the  remote/user  set.  About  20  km  is  the 
current  maximum  range  at  which  the  system  can  operate,  due  to  radio  limitations.  The  real  limit  of 
the  system  has  not  been  established.  According  to  Dr.  Remondi,  under  good  conditions,  it  should  be 
possible  to  initialize  the  system  out  to  40  km  with  the  current  configuration. 

CONCLUSIONS 

Centimeter-level  accuracy  in  three  dimensions  in  real  time  is  a  reality.  The  OTF  system  can  be 
used  on  many  mobile  platforms  and  has  been  tested  on  several  types  under  varying  speeds  and 
dynamic  conditions.  The  tests  have  convinced  us  that  this  technology  will  have  a  great  impact  not 
only  for  surveying  and  navigation,  but  in  other  areas  as  well.  Current  plans  include  exploring  the 
possibilities  in  the  area  of  deformation  monitoring  of  large  engineering  structures  as  well  as  applications 
in  robotics.  This  particular  system  is  robust,  reliable,  and  easy  to  use.  The  tests  have  shown  that  real¬ 
time  modeling  of  tides  using  OTF  is  possible;  this  is  very  important  for  the  dredging  industry.  Although 
this  system  is  referred  to  as  a  prototype  and  will  be  further  developed,  it  has  far  exceeded  the  original 
design  specifications  and  can  be  used  today  as  a  reliable  working  system  for  applications  that  require 
real-time  centimeter  positioning.  It  also  is  a  valuable  engineering  tool  for  those  operations  for  which 
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post-processed  data  will  suffice.  In  reality,  we  have  demonstrated  that  real-time  KOTF  GPS  is  as  easy 
to  provide  as  OGPS  within  the  current  range  limitations. 

PI5CIAIIVIER 

The  mention  of  a  commercial  company  or  product  does  not  constitute  an  endorsement  by  the 
National  Oceanic  and  Atmospheric  Administration  or  the  U.  S.  Army  Corps  of  Engineers.  Use  for 
publicity  or  advertisement  purposes  of  information  from  this  publication  concerning  proprietary 
products  or  the  test  of  such  products  is  not  authorized. 
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